Absence of Translational Ordering In Driven Vortex Lattices 



On 
ON 



O 

On 
(N 



Steven Spencer^ and Henrik Jeldoft Jensen 
Department of Mathematics, Imperial College, 180 Queen's Gate, London, SW7 2BZ, United Kingdom. 

(February 1, 2008) 

By using finite temperature molecular dynamics simulations, we consider the question of the 
existence of dynamical ordering as a two dimensional vortex lattice is driven through a random 
background potential. The shape of the current-voltage(I-V) characteristics are qualitatively the 
same as those seen in experiments on the conventional superconductors 2_H-NbSe2 and Mo77Ge23. 
However, in contrast to previous simulations, we find that the lattice has no topological or trans- 
lational order for any value of the driving but that it does exhibit a high driving phase where long 
range orientational order exists. 

PACS: 74.60.Ge, 64.60.Ht, 05.70.Fh 



o 
o 

Of 



-a 
c 

o 
o 



> 
o 

(N 
O 

On 

-i— > 
c3 



i 

O 
o 



X 



I. INTRODUCTION 



Generic lattices that are driven through a disordered 
environment provide a classic example of an elastic 
medium that is driven through a random quenched back- 
ground potential, and a particular example of these is 
the driven vortex lattice in dirty type II superconduc- 
tors. Other examples are charge density waves in disor- 
dered conductors, fluids in porous media and motion of 
magnetic bubbles through a disordered medium. In the 
former case, the vortex lattice, which is set up when the 
superconducting material is placed in a magnetic field of 
magnitude between its upper and lower critical fields, is 
pinned by randomly positioned defects in the crystalline 
structure of the material. When an external current is 
driven through the material, the motion of the vortices is 
also impeded by the disorder in the material. This results 
in the existence of a depinning transition at zero temper- 
ature, below which the average velocity of the vortices 
(which is proportional to the voltage across the sample) 
is zero, while above it is finite. The depinning transition 
was initially suggested to be a dynamical critical phe- 
nomenon, with diverging length scales at the depinning 
transition, and this was treated for the case of charge 
density waves using mean-field elasticity theory jl| . The 
fact that elasticity theory was used meant that topolog- 
ical defects within the charge density wave (or vortex 
lattice) where neglected. 

Several years ago, the simulations of Jensen and co- 
workers [HH showed that as the vortex lattice was driven 
through a random quenched background potential, it did 
not depin as a single coherently moving elastic medium 
but, rather, the lattice began to break and flow plasti- 
cally. That is to say that there were regions of the lattice 
that were pinned by the random potential and then chan- 
nels of moving vortices would flow around them. Their 
simulations showed that the strength of the random po- 
tential that was required to make the lattice flow plas- 
tically vanished in the thermodynamic limit and, thus, 
elasticity theory would always break down. Even after 



taking this into consideration, elasticity theory has still 
been frequently used to study the depinning transition 
in many systems, particularly the depinning of a d — 1 
dimensional interface in a d-dimensional disordered envi- 
ronment 0| . 

Recently, the focus has turned away from the study 
of the transition itself and turned towards the physics 
of the moving phase. In particular, since the vortex lat- 
tice flows plastically near to the critical driving and thus 
will contain much topological disorder, does the system 
exhibit a dynamical phase transition from a low driving 
phase, which has a high degree of topological disorder, 
to a high driving phase which is topologically ordered? 
This question was first addressed by Shi and Bcrlinsky 
Q and Koshelev and Vinokur Q. Both sets of authors 
performed computer simulations on simple models of a 
two dimensional vortex lattice in a random pinning po- 
tential and both have answered the above question in the 
affirmative. 

Subsequently, Giamarchi and Le Doussal @ have per- 
formed analytic studies of such a system and they pre- 
dicted that the high driving phase of the system is actu- 
ally a moving topologically ordered glass. Subsequent to 
this have been further simulations by various authors, on 
flux line systems [jsj-jTof and on charge density wave sys- 
tems |]ll| , most of which claim that there is a transition 
to a topologically ordered high driving phase. 

The experimental situation is much less clear since it is 
difficult to examine the fine structure of the moving lat- 
tice directly and, thus, indirect methods can only be used. 
Experiments on bulk 2H-NbSe2 Q and Mo77Ge23 films 
Jl3| involve measuring the I-V characteristics and then 
deducing the differential resistance from the data. These 
I-V curves, which are carried out in the so-called peak 
regime 14 1 , bare strong resemblance to early simulations 



of Jensen and their shape has been attributed, by the au- 
thors, to plastic flow. In the experiments of Bhattacharya 
and Higgins , they claim that the dynamical ordering 
in their experiments occurs when the differential resis- 
tance reaches its high driving asymptotic value, while 
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the authors of reference |]13[ claim that this occurs at 
the peak in the differential resistivity. Since the methods 
used to make these claims are indirect, i.e. they are taken 
from the I-V characteristics alone, simulations are a very 
useful tool in studying the physics of the driven system in 
far more detail. The simulations that we present herein 
show that the more correct interpretation is that of ref- 
erence |L2] ], However, even though we show that there is 
a crossover between regimes at the onset of the asymp- 
totic flux flow resistance, the high driving phase is not 
topologically ordered and topological defects are always 
frozen into the lattice structure. 

The rest of the paper is outlined as follows: the model 
and method are presented, followed by the results, dis- 
cussion and finally our conclusions. As per usual. 



II. MODEL AND METHOD 

The model we consider is that of a two dimensional 
array of pancake vortices. The vortices are contained 
within a rectangular box with periodic boundary condi- 
tions. The dimensions of the box are chosen such that its 
size is commensurate with the formation of a perfect tri- 
angular vortex lattice. That is to say that, for a perfect 
vortex lattice of lattice spacing ao containing N v — N x N y 
vortices, the dimensions of the box are L x = N x ao and 
L y = Nyao^/3/2, thus giving an areal density of vortices 
of n v = N v /(L x L y ) = 2/(\/3ao)- ^ n our simulations, we 
set our unit of length by ao = 1 and the system contains 
N v = 34 x 30 = 1020 vortices. The vortices then inter- 
act through a pairwise Gaussian repulsive potential, such 
that the potential energy of the interaction between the 
ith and the jth vortex is 



U vv {ri - rj) = A v exp - 




(1) 



where A v is the strength of the interaction and R v is its 
range. We set the energy scale of the simulations here by 
putting A v = 1. In all the simulations presented herein, 
the vortex-vortex interaction range is set at R v = 0.6, 
which corresponds to a shear modulus of cqq — 0.2695 
and a compression modulus of en = 1.9943 ||. 

The pinning potential in the simulations is generated 
by a set of N p randomly positioned pinning centres with 
positions Ft? . The vortex-pin interaction is then given 
by an attractive Gaussian, such that the potential energy 
of the interaction between the ith vortex and the jth pin 
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-A p exp 
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(2) 



where A p is the strength of the interaction and R p is its 
range. The pinning parameters used in these simulations 



are R p = 0.25 and A p — 0.5. The density of pinning 
centres is n p = N p /N v = 1. This set of parameters en- 
sures that the plastic flow of vortices occurs in the driven 
system. Let us outline the simulation method below. 

From the above model parameters, we can write down 
a Hamiltonian for the system 



H 



N v N p 



r i- r j)+LL^( r i- R r) ( 3 ) 



--1 j=l 



Before we simulate the I-V characteristics of the sys- 
tem, we relax the vortex configuration into the pinning 
potential at zero driving force. Therefore, the initial con- 
ditions for the simulations of the I-V characteristics are 
a relaxed lattice that already has many topological de- 
fects, rather than a perfectly ordered lattice as chosen 
by several other authors [|]](| . We do this relaxation by 
an annealing technique using Newtonian mechanics and 
a rescaling of the velocities after each time step An 
initial temperature of Tj = 1 is chosen and then the sys- 
tem is annealed by slowly reducing the temperature to a 
final value of Tj = 10~ 8 . At the end of this part of the 
simulation, we store the vortex configuration and use it 
as the initial configuration for calculating the dynamical 
properties of the driven vortex lattice. 

From the Hamiltonian above, we can then write down 
the equation of motion governing the over-damped dy- 
namics of the driven system 



di"i dTL „ , . 

at ori 



(4) 



the first term on the right hand side of the above equation 
represents the total force on vortex number i due to its 
interaction with all the other vortices and the pins, the 
second term is the homogeneous driving force applied to 
each vortex F<j r = Fdr*-, which mimics the applied cur- 
rent. These first two terms make up the deterministic 
part of the forces, while the third term on the right hand 
side is a stochastic term, which has the statistical prop- 
erties 

(Xi(t)) = and ( X i(t)Xj(t')) = A6 l3 6(t - t') (5) 

Here, the angular brackets denote averages over the full 
distribution of random variables and the amplitude of the 
two-point correlator can be related to an effective tem- 
perature through A — 2ksTri pj|] . In these simulations 
we set the Boltzmann constant and the friction coeffi- 
cient to Ub = tj = 1 and we use the discrete form for the 
stochastic forces as outlined in ref . M| . 

Firstly, we perform simulations on the vortex sys- 
tem in the absence of pinning and at zero driving force 
(A p ,Fdr = 0). We study the structure of the lattice as 
a function of the temperature to determine the melting 
scenario for the system. We measure the density of mis- 
coordinated vortices, n m i S , derived from constructing the 
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Voronoi Diagram for the vortex array the transla- 
tional order parameter, 



vfr,. 



-c: 



(6) 



and the hexatic order parameter, 



®H = 



1 N " 1 " 3 



,160, 



(7) 



Here, both sets of angular brackets denote time averages, 
G is a reciprocal lattice vector, Zj is the coordination 
number of vortex j and Ojk is the bond angle between 
vortex j and its kth neighbour relative to the fixed x- 
direction. When we calculate these quantities for the 
driven system, we always evaluate the quantities in the 
co-moving frame. 

For a two dimensional lattice in the thermodynamic 
limit, translational order only exists at T = 0. For any fi- 
nite temperature, the long wavelength modes destroy the 
translational order and — 0. However, even though 
long range order is absent, the translational correlation 
function has a power law form rather than the exponen- 
tial form of an isotropic liquid. Thus, the system is said 
to have quasi-long range order. In continuum elasticity 
theory, one can write the mean squared displacements of 
the vortices as 



(u 2 ) th 



2ir 2 



dk 

c 66 k 2 . 



(8) 



where the limit en S> cqq has been assumed. Thus, one 
can see that the divergence at k = causes unbounded 
displacements which gives the loss of translational order. 
One can prevent this divergence by introducing a lower 
cut-off in the integration brought about by considering a 
finite sized system of linear extent L. The Brillouin Zone 
then spans the A-vectors 7r/i < k < n/aQ and one can 
perform the integration to give 



/ 2\ _ k B T L 

(u )th ~ m — 

C667T ao 



(9) 



i.e. bounded displacements. One can then make an es- 
timate of a melting temperature by applying the Lindc- 
mann criterion, i.e. when the displacements are u ~ clclq 
(where cl is the Lindemann number), the lattice has 
melted. This gives a melting temperature, T m , of 



k B T, 



B -L m 



c 6e TT(c L a y 
ln# 



(10) 



which vanishes as L — > oo. One can invert the above 
equation to give L or {T). Here, L or is the linear size of 
a region (in an infinite system) over which translational 



order is maintained, and one can even think of the vol- 
ume l? or as being analogous to the correlated volume in 
collective pinning theory JlT]] , with the size of the volume 
diverging as T — > 0. However, the above equation for the 
"melting" temperature (written in " " because it is re- 
ally just a cross-over where L or becomes of the order of 
the system size) has only a very weak inverse logarithmic 
dependence on L and thus we can expect that for macro- 
scopic, but finite systems, a low temperature "phase" 
exists where translational long range order exists (again 
written in " " since it is not a true thermodynamic phase 
transition). 

On increasing the temperature, this two dimensional 
solid then undergoes a melting into a liquid phase. In 
the Kosterlitz-Thouless type scenario for two dimensional 
melting, the two dimensional solid with quasi-long range 
order melts at a temperature T^ T = 0.62c66Oq/ (An) p8[ , 
into a hexatic liquid (for the system that we consider, 
this is T^ T w 0.0133). This melting is brought about by 
the unbinding of dislocation pairs and the hexatic liquid 
is characterised by short range translational order (with 
= as before) and quasi-long range orientational or- 
der, i.e. ^ h = with power law hexatic correlations. 
At a higher temperature still, the unbound dislocations 
undergo a further unbinding by splitting into unbound 
disclination pairs. Here, both the translational and ori- 
entational order are short ranged and an isotropic liquid 
is formed. 

We then simulate the temperature dependence of 
n m is , and m order to compare with the above 
scenario. Figure [l] shows the three quantities along with 
the structure factor S(k) for two different temperatures. 
One can clearly see that there is a finite temperature, T m , 
at which both the translational and hexatic order param- 
eters vanish. This also coincides with the rapid increase 
in the density of miscoordinated vortices in the system. 
The structure factors S(k) are shown for temperatures 
just above and below T m . The system clearly undergoes 
a transition from a solid, signified by the set of sharp 
delta-function peaks, to an isotropic liquid, signified by 
the single k — peak and concentric rings. The peaks in 
the solid with quasi-long range order diverge algebraically 
while in the hexatic liquid, 5(k) should display a set of 
concentric rings with six- fold symmetry |l^,^] . In these 
simulations, it seems that there is a true solid phase for 
low temperatures, i.e. translational and orientational or- 
der exist throughout the considered system (L or > L), 
which then melts into an isotropic liquid at a finite tem- 
perature T m . Thus, no solid with quasi- long range order 
or hexatic liquid are observed. 

From these simulations, the melting temperature is 
then T m w 0.011, where we have used the structure fac- 
tors and order parameters as signitures of the high and 
low temperature phases. 
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FIG. 1. The pure, undriven system - the density of mis- 
coordinated vortices and the translational and hexatic order 
parameters. Also shown are the structure factors, 5(k), slight 
below and slightly above the melting transition. 



From these simulations, we can choose an upper bound 
for the temperature at which the driven simulations are 
to be run, namely for temperatures below T m . This is 
the obvious upper bound since, when investigating the 
dynamical ordering of the vortex lattice as it is driven 
through the pinning potential, there can clearly never be 
any ordering for temperatures above T m . 

We now return to the system that contains pinning 
in order to address the central question of this article. 
When driving the lattice through the random potential, 
we solve the set of equations |] numerically for a partic- 
ular value of the driving force. The numerical time step 
for the simulations is chosen as that which gives good 
energy conservation in the case of Newtonian dynamics 
(without velocity rescaling). Obviously energy will not 
be totally conserved in this case due to the numerical ap- 
proximation. Instead, it has fluctuations around a con- 
stant average value. We chose a time step of St = 0.01, 
which gives a energy fluctuations of about 0.01% of the 
average. The threshold driving force at zero temperature 
for the parameters stated is F c w 0.545. We begin the 
simulations at low driving Fd r = 0.75 and from the disor- 
dered configuration obtained after the initial relaxation 
procedure. The equations of motion are numerically it- 
erated until the system has had time to respond to the 
change in driving and a steady state has been reached. 
For this condition to be fulfilled, we find that it is more 
than sufficient to solve the equations of motion until the 
centre of mass position of the initial vortex configuration 
has moved 20 lattice spacings. After this initial period, 
time averages are performed on the system. Time aver- 
ages are performed over a minimum of 10000 time steps 
and until the averages have become stable to within a 
factor of ±10 -5 . This always corresponds to centre of 
mass motion of at least 2 lattice spacings, and is usually 
much more. 



III. RESULTS AND DISCUSSION 

Here we present the results of our simulations for 
the set of system parameters listed above. The simu- 
lations where performed at temperatures of T = 0.0001 
and T = 0.001, both below the melting temperature of 
T m » 0.011. 

The four graphs in Figure |^ show the differential re- 
sistance Rd — d{v) /dFdrj t ne fraction of vortices with 
average velocity (v) — 0, N stop /N v , the fraction of mis- 
coordinations in the lattice, N m i S /N v , and the average 
total pinning force, F p , all as a function of driving force, 

Fdr- 
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FIG. 2. a) The differential resistance, b) The fraction of trapped vortices, c) The fraction of miscoordinated vortices, d) The 
average total pinning force. Results for two different temperatures, T — 0.0001 and T = 0.001, are shown. 



The differential resistance curve is strongly resemblant 
of that measured in experimental systems such as 2H- 
NbSe2 12 and those in Mo77Ge23 thin films |l3|], where 



the I—V characteristics have been measured in the peak 
regime JlJ]. The peak in this curve, or the fact that 
R™ ax > Rd{Fdr — * oo) is indicative of the existence of the 
prominent shoulder feature in the I—V curve. If one ex- 
amines the vortex flow patterns at driving forces around 
the peak in Rd, one indeed observes plastic vortex flow in 
channels, flowing around regions of pinned vortices. This 
is in agreement with the comparison between experimen- 
tal curves of this form and this type of vortex flow jl2) . 
Below we present results that show that the only fea- 
ture of the dynamics of the system that can be extracted 
from this particular feature of the I—V curve is that plas- 
tic flow is occurring in the vortex dynamics. One cannot 
determine whether or not dynamical ordering takes place 
in the system from this feature alone. 

From figure ||a, one can see that the peak in Rd oc- 



curs at approximately Fdr ~ 1-8. From figure |^b, one 
can also see that the number of trapped vortices van- 
ishes only above a driving of Fd r ~ 2.8. From this figure 
we can easily conclude that, for this particular system, 
plastic flow definitely occurs for driving forces Fd r < 2.8. 
It should also be noted that this is a lower bound on 
the threshold for plastic flow since the true definition of 
plastic flow is if the lattice is continually being plastically 
deformed as it flows, i.e. the vortices are always flowing 
plastically if some of the vortices are constantly changing 
their neighbours. We can see where the true threshold 
for plastic flow is by considering the distributions of the 
time-averaged individual vortex velocities, which are pre- 
sented below. 

Figure shows the density of defects as a function 
of driving force, with the main features being the strong 
peak at a driving force slight greater than threshold and 
the asymptotic finite value of the density for large driv- 
ing forces. This shows clearly an absence of a dynam- 
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ical transition to a topologically ordered state and this 
absence is a result of hysteretic effects in the system. 
Depending on the initial conditions used in the simula- 
tions, one can obtain the results presented above or one 
can obtain the results of various other numerical works, 
namely that two distinct phases exist - a topologically 
ordered high driving phase and a disordered low driving 
phase which may or maybe not separated by a dynam- 
ical phase transition. Now it remains for us to decide 
which initial conditions arc more viable. In the simula- 
tions of Shi and Berlinsky || and Koshelev and Vinokur 
they began their simulations from an ordered lat- 
tice at large driving force. On reducing the driving, they 
then see a sharp increase in the density of defects within 
the vortex array. Indeed, we can observe this sort of be- 
haviour using the same method. However, the results we 
present above correspond to simulations with different 
initial conditions. We start with a lattice configuration 
that is initially annealed into the pinning potential and 
contains a high degree of topological defects (approxi- 
mately 50% of the vortices are miscoordinated) . The 
simulations then start at a driving force slightly greater 
than the threshold driving force which is then increased 
until the dynamics reach their asymptotic behaviour. We 
suggest that our initial conditions are more realistic and 
correspond to something along the lines of field-cooling 
a sample and then measuring the I-V curves by increas- 
ing the external current. Thus, in these simulations we 
see that, rather than moving into a topologically ordered 
state on increasing the driving, the system moves into a 
state where the number of defects in the lattice become 
constant. It is worth noting that in ref. |L3]], the authors 
point out that a levelling off of the differential resistance 
at high currents is concurrent with a constant density of 
defects. This is indeed the case as shown by our simula- 
tions, with both quantities reaching their asymptotic val- 
ues at around the same value of driving. It is also worth 
noting that other recent simulations [|10| also observe a 
non-vanishing density of defects on increased driving. 

The final graph, figure ||d, shows the average total 
pinning force as a function of Fd r - Since the physical 
quantities presented in figures ||a-c are all a direct conse- 
quence of the existence of the pinning force, the pinning 
force's apparent insensitivity to temperature explains the 
insensitivity to temperature of the other quantities. We 
think nobody should expect that, if the pinning force 
was temperature independent, strong temperature de- 
pendence should be observed in the other quantities. For 
driving forces greater than Fd r ~ 2.0, the pinning force 
is independent of temperature and, in this region, one 
can fit the behaviour to (F p ) ~ -F^ Q , where the expo- 
nent is a rs 9/8. This compares well with an exponent 
of unity obtained from perturbation theory [Tc| ]. The 
fact that (F p ) has a power law dependence may have 
consequences in the thermodynamic limit, which will be 
discussed later. 
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FIG. 3. The same as in figure 1, but for the driven system. 
The structure factors are shown for phases with no hexatic 
order and for hexatic order only. 



Thus, we have seen that no topological ordering takes 
place. This does not immediately imply there is no or- 
der whatsoever within the lattice as we know from the 
Kosterlitz-Thouless scenario. Figure shows the den- 
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sity of miscoordinations, the two order parameters and 
the structure factors for the vortex system at a temper- 
ature of T = 0.0001 and for various driving forces. The 
translational order parameter, represented by the filled 
circles, remains at ~ for all driving forces, but 
from the hexatic order parameter, we can see that, for 
high driving forces, the system has some orientational 
order. This can also be seen from the two structure fac- 
tors shown. For Fd r — 1.7, S(k) just consists of a single 
peak at k = with no other structure. For the large 
driving forces, other peaks have appeared. The peaks 
are anisotropic corresponding to the anisotropic nature 
of the order of the lattice. We are unfortunately unable 
to extract the functional form of the decay of the peaks. 
Thus, there is a dynamical ordering that takes place, but 
it is not to a topologically ordered state. Finally we no- 
tice that the stucture factor for Fdr — 7.0 is representa- 
tive of the structure factors found at any stronger driving 
force. The peaks do not appear to become more narrow 
with increased driving. 

Figure ^| shows a set of graphs showing the distribution 
of the time- averaged individual vortex velocities, with 
each graph corresponding to a different value of the driv- 
ing force. For each of the N v vortices, we calculate the 
following 
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i.e. the angular brackets denote a time average and the 
subscript i denotes the fact that a spatial average has not 
been performed. The figures in || then show the distribu- 
tions of these quantities. Looking at the figures, one sees 
precisely what one would expect for plastic flow occur- 
ring in the system. Slightly above threshold (Fdr = 0.9) 
the distribution consists of a strong peak at (v) = 0, with 
some small but finite velocities having lower weights, i.e. 
most of the vortices are trapped for this value of the 
driving, with the few moving ones contributing towards 
the average centre of mass velocity of the whole vortex 
system. The following graphs for d riving Fdr ^ 2.8 are 
qualitatively similar, with the peak at (v) = decreasing 
and a second peak at finite velocities developing as the 
driving increases. This continues until Fd r ~ 2.8, where 
the peak at (v) = vanishes and all the vortices become 
depinned. This is consistent with the result in figure ||b. 
The figures corresponding to Fd r > 2.8 thus consist of 
a broad distribution with a single maximum. The posi- 
tion and the height of the maximum of the distribution 
both increase correspondingly with the driving force, and 
the width of the distribution, which we call A(i>), clearly 
decreases on increased driving. The fact that these dis- 
tributions are of finite width (A(v) ^ 0) is the singular 
most important fact in determining whether or not the 
vortex array is flowing plastically and is therefore of im- 
portance when addressing the question of whereabouts 
dynamical ordering occurs. 
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FIG. 4. The distribution of time-averaged velocities 
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For the vortex array to undergo a dynamical transi- 
tion to a topologically ordered state at some value of the 
driving force, all free topological defects within the lat- 
tice must vanish and the lattice must also clearly return 
to the regime where only elastic deformations exist. If 
plastic flow exists then topological defects will be con- 
stantly appearing and disappearing. Thus, as we have 
seen from figure ||c that topological defects are always 
present, we need to tie in all the results presented above 
to develop an overall picture of the dynamics. From fig- 
ure |f3 we have stated that, for Fd r < 2.8, some regions 
of the vortex lattice are permanently trapped while the 
other vortices flow aroung them. Thus for driving forces 
in the range F c < Fd r < 2.8, plastic flow occurs with 
some of the vortices being permanently pinned. This is 
in contrast to the work of ref. |l(J which states that some 
of the vortices only remain pinned for a finite length of 
time. If one was to plot the instantaneous distribution of 
vortex velocities, one would see the double peaked struc- 
ture, similar to ours, but on plotting the distribution of 
time-averaged velocities, the double peaked form may dis- 
appear or simply the = peak must move to a finite 
velocity. The fact that the distributions presented above 
are time-averaged and show a double peaked structure 
with one peak at (uj) =0 can only lead to the conclusion 
that these simulations show that the channel structures 
that are formed, along which the plastic flow occurs, are 
'frozen' in. This may, of course, not necessarily be a 
unique structure dependent only on the value of the driv- 
ing force, but probably will be dependent on the history 
of the system. 

Therefore, we know that, for Fd r > 2.8, all of the vor- 
tices are depinned, i.e. (Wj) > Vi. However, this does 
not imply that plastic flow has ceased. Here we return 
to the distributions of figure ||. For all plastic flow in the 
vortex dynamics to be truly eradicated, the infinite time- 
averaged velocity distribution must be a delta function, 
i.e. all velocities are the same and the system must be in 
the elastic regime. Obviously this is not the case for the 
time-averages in the simulations. For finite time-averages 
the width of the velocity distribution can be finite and 
still no plastic flow occurs over the time span. Thus, if 
the average is performed over a time T and the width of 
the velocity distribution is A(v), then one can say that 
the corresponding width of the distribution of distances 
moved during the time span of the average has a width 
of Ax w A(v)T. We can then say that, for no plas- 
tic flow to occur in the simulation, Ax < oq/2. Thus, 
say for a time average performed over 10000 steps (re- 
member St = 0.01), plastic flow will always occur in this 
time window if the width of the velocity distribution is 
A(v) > A(v) p i = 0.005. For the simulations presented 
herein, the time-averages are always performed over a 
time T > 10000<ft = 100 and, thus, A(v) p i = 0.005 is 
an upper bound for our simulations, i.e. any distribution 
width taken from these simulations that is greater than 



this corresponds to plastic flow. Thus we are in a posi- 
tion to tell when then vortex system has returned to the 
elastic regime. 

From figure [|, one can see that the width of all the 
velocity distributions is greater than this threshold for 
plastic flow. In fact, only as the driving force Fd r ap- 
proaches a numerical value of w 10 does the width A(v) 
become less than the threshold. Thus, to bring together 
all the results above, we can present the following picture. 

The dynamics above threshold can be split into two 
'phases', with one of these phases being split into two 
'sub-phases'. Figure [| shows the scenario that we present 
as a schematic diagram. The first phase is for driving 
forces in the range F c < F dr < F 2 (F c » 0.545, F 2 » 10), 
where the vortices flow plastically and the phase can be 
defined through a non-zero width of the time-averaged 
individual vortex velocity distribution (presented in fig- 
ure ||), A(v) > 0. This phase is then split into two sub- 
phases. These are separated by a threshold F\ « 2.8, 
which is the second depinning force. The first depinning 
force is simply that which separates phases of zero and 
non-zero centre of mass velocity, i.e. F c . F\ signifies 
the point where all of the vortices become depinned. For 
F c < Fd r < F\, some vortices are moving but some are 
permanently pinned ((v)i = 0) and thus N stop > 0. For 
F\ < Fdr < F 2 , all the vortices have finite time-averaged 
velocities but the width of the velocity distribution is 
finite, signifying plastic flow. 
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FLOW 


FLOW 
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N sto ^>0 j N sto =0 











FIG. 5. The dynamical 'phases' of the driven vortex sys- 
tem. For the simulations presented herein, F c ~ 0.545, 
F\ ~ 2.8 and Fa ~ 10. See the text for a more comprehensive 
explanation. 

The second phase is for Fdr > F 2 and is defined 
through A(v) = 0. Here, all plastic flow has ceased and 
the vortex array moves as a single coherent structure. 
That is to say that all the vortices keep the same neigh- 
bours throughout the simulation. Here we use the phrase 
'single coherent structure' as opposed to saying that the 
system has returned to the elastic regime in order to avoid 
any association the reader might make between an elastic 
regime and a topologically ordered structure. 

Given that these are the different regimes of dynam- 
ics in these simulations, we now compare how the time- 



8 



averaged quantities vary in each of these regimes. The 
quantities of interest are the density of miscoordinations, 
the two order parameters and the differential resistance. 
Firstly, looking at figure ||a, one can see that the only fea- 
ture of the differential resistance that correlates with the 
regimes described above is that the onset of the asymp- 
totic flux flow value of the resistance coincides with the 
third regime (starting at Fd r ~ 10), i.e. where the plas- 
tic flow in the system ceases. A similar thing can be said 
about the density of defects. The flattening off of the de- 
fect density also coincides with the onset of this regime. 
Thus, our main statement of the work is the following: 

There is a lack of a dynamical transition to a topologi- 
cally ordered state in two dimensional vortex lattices with 
strong pinning. However, there are two possible dynam- 
ical regimes that can be distinguished. The first is a low 
driving regime, where the vortices flow plastically defined 
through A(v) > (pinned regions may or may not exist). 
Here the differential resistance and the density of defects 
vary strongly with the driving force. The second regime 
is the high driving regime. Here, the vortex array moves 
as a single coherent structure (A(u) = 0), i.e. all plastic 
flow has ceased, and the vortex dynamics have reached 
their asymptotic behaviour, i.e. the differential resistance 
and the density of defects are independent of the driv- 
ing force. Since all plastic flow has ceased in this phase, 
the defects that remain in the lattice are frozen in. This 
is also reflected in the fact that the translational order 
parameter is zero for all driving forces. Thus, the high 
driving phase is not a topologically ordered phase, but 
one where the vortex lattice is a single coherently moving 
topological glass. However, even though we have shown 
that no topological or translational order exist at any 
driving force, for large enough driving forces the system 
exhibits long range orientational order, given by 0. 

IV. CONCLUSIONS 

Let us now consider the discrepancy between the sim- 
ulations in this article and some recent simulations by 
other authors |i|||||-[l0) ■ The majority of these simula- 
tions see quite a sharp dynamical transition to a topo- 
logically ordered state at some value of the driving force, 
which is in complete contrast with the simulations herein. 
We suggest that this is either an artifact of the weak 
pinning strengths that these simulations have used or 
an artifact of the initial conditions, as mentioned pre- 
viously (the topological ordering observed by some au- 
thors is truly an artifact of their initial conditions since 
they are rather unrealistic initial conditions). Here, like 
most other authors, we are guilty of making a distinc- 
tion between weak and strong pinning. In the thermo- 
dynamic limit this distinction becomes meaningless. In 
the early simulations of Jensen and co-workers they 
showed that the pinning strength required to produce 



plastic flow in the system vanishes logarithmically as the 
thermodynamic limit is approached and, thus, the sys- 
tem is always in the plastic flow regime . We suggest 
that the same effect could happen here and, thus, the 
dynamical ordering seen in other simulations is simply a 
finite size effect. 

We have also stated here that, on increased driving, our 
simulations return to a regime where plastic flow ceases. 
Admittedly, this could well be a finite time effect since we 
have stated previously that, in the limit of infinite times, 
the distribution of time-averaged velocities must be a 5- 
function for all plastic flow to cease. Thus, only when 
the width of the distribution does vanish will plastic flow 
truly cease. It is also worth noting that this may never 
occur. We have shown that the average pinning force 
never vanishes, i.e. it behaves as (F p ) ~ 1/Fd r - This 
shows that the asymptotic flux flow resistivity is never 
actually reached in this model and, for infinite systems 
and times, any finite pinning force is capable of producing 
plastic flow. 

Thus, we suggest the following. The simulations in 
this article have still not completely resolved all the ques- 
tions as to dynamical ordering in vortex lattices. What 
is needed is a detailed finite size study of the dynamical 
regimes and topological phases of the lattice so that a 
coherent picture of the dynamical behaviour for the infi- 
nite system can be evaluated. We have shown that, for 
our system sizes, topological and translational order do 
not exist and we believe that considering larger systems 
will not change this 0. However, we have observed a 
low driving phase without orientational order and a high 
driving phase with orientational order. From a purely 
theoretical point of view, it would be interesting to see 
whether or not these phases persist in the thermody- 
namic limit and whether or not they are separated by 
a true dynamical transition. We have suggested, above, 
that plastic flow may persist for all driving forces in the 
infinite lattice. From the figures presented herein, we can 
also see that the onset of hexatic order overlaps with the 
region where plastic flow occurs. Thus, we can suggest 
that plastic flow does not destroy the hexatic order and 
that this order could easily persist in infinite systems. 
Once the physics of the infinite system is realised, the 
discrepancies between our simulations and many others 
could probably be easily resolved. 
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